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Abstract 

The one-dimensional t-J model is investigated by the variational Monte Carlo 
method. A variational wave function based on the Bethe ansatz solution is 
newly proposed, where the spin-charge separation is realized, and a long- 
range correlation factor of Jastrow-type is included. In most regions of the 
phase diagram, this wave function provides an excellent description of the 
ground-state properties characterized as a Tomonaga-Luttinger liquid; Both 
of the amplitude and exponent of correlation functions are correctly repro- 
duced. For the spin- gap phase, another trial state of correlated singlet pairs 
with a Jastrow factor is introduced. This wave function shows generalized 
Luther-Emery liquid behavior, exhibiting enhanced superconducting correla- 
tions and exponential decay of the spin correlation function. Using these two 
variational wave functions, the whole phase diagram is determined. In addi- 
tion, relations between the correlation exponent and variational parameters 



in the trial functions are derived. 
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I. INTRODUCTION 



The anomalous properties found in high-Tc superconducting copper oxides have led 
to a renewal of interest in strongly correlated electron systems in low dimensions. Among 
various candidates, the t-J model has attracted considerable attention as a model to describe 
the cuprate superconductors. |^,^ 

For the one-dimensional (ID) t-J model, much progress has been achieved using ana- 
lytical and numerical techniques. |Q It has been found [^-0] that the three main regions 
can be distinguished in the phase diagram defined by the electron density and the ratio of 
spin exchange interaction to hopping amplitude, J/t. First, a Tomonaga-Luttinger liquid 
(TLL) holds for small J/t, which is characterized by power-law decay of correlation 
functions. It has been clarified that the separation of spin and charge degrees of freedom is 
playing an essential role in this region. [0 Second, phase separation takes place for large 
J/t, where the system is separated into electron-rich and hole-rich phases. Third, there is a 
region with a gap in the spin excitation spectrum for J /t > 2 and at small electron densities. 



On the other hand, in the 2D t-J model, although some aspects are obtained so far, |11 



many problems are left unresolved. Particularly, the crucial question is whether the features 
realized in ID system, like the charge-spin separation and/or TLL, take place also in 2D 
system or not. To obtain an unified and consistent understanding of the 2D t-J model, 
further progress of research is needed. 

The variational Monte Carlo (VMC) method is one of the most powerful and transparent 
approaches to investigate strongly correlated electron systems. It provides a deeper 
insight because of its explicit form of the wave function. It is very important to construct 
a better trial function in the framework of the VMC technique. One way to obtain further 
insight into the wave function in the 2D t-J model is to extend the wave function realized in 
ID system. For this purpose, examining trial wave functions for ID system in detail gives 
us useful references in the pursuit of the 2D t-J model. We shall study variational wave 
functions in the ID t-J model, keeping a possibility of extending to the 2D system in mind. 
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So far, various kinds of variational functions have been proposed for strongly correlated 
electron systems. p!^^l9[ The Gutzwiller wave function ||20[ has been extensively studied for 



its simplicity, and shown to be a good trial function for the supersymmetric (J/t = 2) ID 
t-J model. This wave function was improved for other values of J/t by introducing a 
conventional Jastrow-correlation factor, but the expected TLL behavior was not recovered. 
17| , p!9|| Recently, Hellberg and Mele have introduced a simple trial wave function of Jastrow- 



type. fl^ It takes into account an effect of long-range correlations, and shows successfully 
the anomalous power-law behavior in correlation functions. This wave function has been 
extended to the 2D t-J model to discuss the TLL instability. |^ However, the properties 



of the exact ground-state are not wholly reproduced by this wave function, especially in the 



small J/t region. ]T9| 



In this paper we introduce another type of variational wave functions. Two kinds of trial 
functions are newly introduced. First, we consider a variational wave function based on the 
Bethe ansatz solution. In the limit of J/t ^ 0, the charge and spin degrees of freedom are 
completely separated, and the ground-state wave function obtained from the Bethe ansatz 
can be written as a product of the two contributions. For finite values of J / 1, although 



the charge and spin are separated, they interact strongly. To take into account this effect, 
we introduce a Jastrow-type correlation factor into the Bethe ansatz solution for J/t —* 0. 
It is shown that this wave function has an advantage of providing an excellent description 
of the ground-state properties in most regions of the phase diagram; Both of the magnitude 
and exponent of correlation functions are correctly reproduced, and a quantitative discussion 
can be made. 

Next, we consider a trial function for small electron densities. For the spin-gap phase, 
Chen and Lee have proposed a variational function of a gas of noninteracting bound singlet 
pairs. [0 This wave function corresponds to a Luther-Emery state p2|,p| with infinite cor- 
relation exponent. More accurate trial function can be generated by correlating the singlet 
pairs with a Jastrow factor. This is just our trial wave function for the spin-gap phase 
introduced in this paper. This wave function shows generalized Luther-Emery liquid be- 



havior, exhibiting enhanced superconducting correlations and exponential decay of the spin 
correlation function. 

Comparing energies of the trial function based on the Bethe ansatz solution and the 
generalized Luther-Emery state, the entire phase diagram is determined. Evaluating the 
correlation exponents by a finite size scaling analysis, the relations between the correlation 
exponent and variational parameters in these trial functions are derived. 

This paper is organized as follows. In the next section our trial functions are introduced. 
Section III provides the results of physical quantities by the VMC calculations. Energies 
and various correlation functions are compared with the exact calculations in Sec. Ill A. 
In Sec. Ill B, the correlation exponents are evaluated from the finite size scaling to discuss 
the long-range behavior of correlation functions. The phase diagram of the ID t-J model 
determined by our wave functions is shown in Sec. Ill C. Section VI is devoted to a summary 
and discussions on related problems. 



II. TRIAL WAVE FUNCTION 

The t-J model is defined by the Hamiltonian, 

H = -tY: (4c,. + H.c.) + J^(S, ■ S, - ln,n,), (1) 

(ijV {ij) 

where cl^. = c|o.(l — nj cj^. being the creation operator for an electron with spin projection 
a at lattice site i, and rii = J2a "-jo- = J2a cJo-Cjo-- Thus creates an electron only on an 
empty site, avoiding double occupancy. The spin operator associated with site i is defined as 
Sj = I J2a,f3 Cia^a,i3Cif3, whcrc <T = {ax, CTy, az) is a vector of Pauh matrices. The summations 
in Eq. (|1]) are taken over nearest neighboring pairs. This model reduces to the U = oo 
Hubbard model in the limit J/t — 0. 

For highly correlated electron systems, Gutzwiller-Jastrow-type wave functions with two- 
body correlation factor are fairly common. HHH The Gutzwiller wave function, ||20| which 



is a prototype of the trial function of this type, is often used as a starting trial function for 
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its simplicity. It is defined as 

IV^g) = Pd |0f) = n (1 - ^^T^^i) I0f) , (2) 

i 

where 0f is a simple Fermi sea and Pd is the operator projecting out the double occupancy. 
This wave function is essentially a Fermi-liquid state, having a discontinuity in momentum 



distribution at k = kp- Thus, the expected TLL behavior was not recovered. Hellberg 



and Mele have introduced a variational state with a long-range correlation: [ITB 



I^hm) = n n [1 - (1 - \d^Jn n^.n,,'] \iPg) , (3) 

dij = sin [7irij/Ns] , (4) 

where rij = jr, — rj\ is the distance between the i-th and j-th sites, and A^^ the number 
of sites. When z/ = 0, ^ipnu reduces to ipG- It has been shown that ipnu exhibits the 



characteristic behavior of TLL. However, the correlation exponent estimated with this 
wave function does not coincide with the exact value for small J/t. |[T^ This disagreement 
becomes apparent when the global features of various correlation functions are compared 



with the exact ones. 19 



An important feature of the TLL is the separation of spin and charge degrees of freedom 
in the low-energy excitations. In the limit of J/t ^ 0, Ogata and Shiba have shown that 
the ground-state wave function obtained from the Bethe ansatz has a simple form due to 



the complete decoupling of charge and spin degrees of freedom. |10| It can be written as a 
product of a Slater determinant of spinless fermions describing the charge degrees of freedom 
and the spin wave function of the squeezed Heisenberg model in which all empty sites are 
omitted. The ground-state wave function in the limit of J/t is expressed as 

ipo{xi, . . . , xnJ = X{xi, . . . , XNjY{yi, ....Vm), (5) 

where 

. . . , xatJ = det [exp(igjXj)] , (6) 
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{xj} are the positions of A^'e electrons, and {yj} are the coordinates of M up spins with 
vacant sites omitted. X{xi, . . . , xjvj is the wave function for noninteracting spinless fermions 
with momenta {qi}, and Y{yi, . . . ,yM) is the ground-state wave function of the S = ^ 
antiferromagnetic Heisenberg modeL 

For finite values of J/ 1, the charge and spin degrees of freedom are no longer completely 
separated, and the charge-spin coupling occurs. Thus we introduce a Jastrow-type correla- 
tion factor in our trial wave function to mix the charge and spin degrees of freedom. We 
shall study the following variational state for the ID t-J model: 

I^ba) = Fj \XY) , (7) 

where the amplitude of |X) = Ylf" cj^ |0) is given by Eq. (||), ensuring the absence of double 
occupancy. The long-range correlation factor of Jastrow-type in Eq. is defined as 

= n n {1 - [1 - Vinr, aa')] n,^nj^,}, (8) 
and the form of function r] is assumed to be 

I Idijl-"^ , if cT = d' 
riirij] aa') = < , (9) 

[ if a 7^ a' 

where dij is given by Eq. (Q). The Jastrow factor Fj modulates the Bethe ansatz wave 
function by the distance between all pairs of particles. Positive value of z/i (z/2) induces a 
repulsive correlation between particles with the same spins (opposite spins), while the neg- 
ative values provide an attractive correlation. When I'l = 1^2, Fj reduces to the correlation 
factor studied by Hellberg and Mele. 

The wave function of the spin part in Eq. is approximated as a trial function of 
Jastrow-Marshall-type, 



F(i/i,...,yM) = (-i)^^^'^ni^^.r% (10) 

where Sij = sm[n{yi — yj)/Ne] for a system of electrons and L {yi} is the number of 
up spins in one sublattice contained in the spin configuration {yi, . . . ,yM)- With this trial 



function Y{yi, . . . jI/m), we have calculated the ground-state energies of the ID antiferro- 
magnetic Heisenberg model, i^Heis = J^{ij) • Sj with N^. < 70 by the VMC technique, 
and estimated the energy in the thermodynamic limit from finite size scaling with a formula 
E/Ne = Eoo + C/N^. The minimum energy is realized for Ug ^ 2, and the resultant energy 
is Eoo = (—0.4421 ± 0.0001) J, which is quite close to the exact value by the Bethe ansatz, 
H EBA/Ne = -(In 2 - \)J = -0.443147. . .J. The difference is only 0.24 %. Therefore, 



Y{yi, . . . jUm) well reproduces the true ground-state wave function of the ID Heisenberg 
model. 

As a result, we have three variational parameters in our trial state (^, i.e., z/i, 1/2, and 
Ug. In most regions of the phase diagram, this wave function ip-^A successfully reproduces 
the exact ground state of the ID t-J model as shown in Sec. III. 

However, a Luther-Emery liquid behavior, exhibiting a gap in the spin excitation spec- 
trum and enhanced superconducting correlations, is found for J > 2t and at small densities. 
I^J^ The true ground state for this region lies out of the variational subspace spanned by 
ipBA- To represent the spin-gap phase better, we introduce another trial state as follows. 

I^rvb) = FjP, J2 n h''"'"~'Kjn] |0) , (11) 

where [i,j] = {cl^Cji~cliCj^) is a singlet pair in a given configuration {injn}, and Pd projects 
out the double occupancy. In Eq. (p^D , controls a weight for a singlet bond as a function 

of its length. The function t] in the Jastrow factor is taken to be rjijij^aa') = \dij\^, i.e., 
Fj is assumed to be spin-independent. Two variational parameters, A and h, are contained 
in the trial function V'rvb- This is a natural generalization of the wave function of a gas 
of noninteracting bound singlet pairs proposed by Chen and Lee, which corresponds to 
a Luther- Emery state with infinite correlation exponent. Correlating the singlet pairs with 
the Jastrow factor Fj, ^/'rvb can be expected to exhibit generalized Luther-Emery behavior. 
It is also a particular form of the resonating valence bond (RVB) state. In fact, ipKVB can 
be rewritten as 
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I^rvb) = FjPa 



^ cos k — h II 



Ne/2 



|0), (12) 



which exphcitly shows a singlet hquid picture of the RVB state. 

III. SIMULATION AND RESULTS 

In this section, we present the resuhs of the VMC calculations for various values of J/t 
and Tie = Ne/Ns with A^e and A^^ being the number of electrons and sites, respectively, and 
make comparisons with those of exact calculations and other trial functions. We consider the 
ID t-J model with up to 300 sites under the periodic boundary condition with Nf,/2 = odd. 

Variational parameters in Eqs. (|^) and (|TT|) are optimized using a conjugate- gradient 
method combined with the fixed sampling in the VMC calculations. Technical details of the 



optimization procedure were described in Ref. [^, and some practical improvements are 
made to achieve the convergence rapid enough to handle multiparameter optimization: a 
quasi- Newton algorithm is employed instead of Powell's optimization algorithm in Ref. ||25| , 



and the gradient is evaluated by the numerical differentiation. Once the fully optimized 
wave function is obtained, we use it in evaluating the physical quantities with another VMC 
run in order to examine the properties of the ID t-J model in detail. Calculated quantities 
are the total energy per site, the momentum distribution function, 

M^) = ^Ee^'^""''^^(4c,.)> (13) 

and the equal-time correlation functions, where (■ ■ ■) indicates the expectation value for a 
given trial function. The spin- and charge- correlation functions in Fourier-transformed form 
are defined as 

Sik) = ^Y.^'"^'''~''^{StS^), (14) 

s ij 

cik) = ^ E ^^'^'■'"'■^^ [M - (^^) (^.)] ' (15) 
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respectively. The singlet pairing correlation function is defined as 

m = ^T.^''^''-"'^ (M^^) , (16) 

where Aj is the annihilation operator of a nearest neighboring electron singlet pair, 

= -^(CitQ+n - Ci|Q+it). (17) 

We collect typically 30,000 samples to take averages of the energy for the optimization of 
variational parameters, and 100,000 - 200,000 samples for the evaluations of the expectation 
values of observables. 



A. Quarter-filled case 

First we compare the properties of ipBA with those of ^ipHM for the quarter-filled case: 
Ue = |. At this electron density, the spin-gap state is absent. In addition, when J/t ^ 0, 
the exact results of correlation functions have been obtained ||T0| , p!7|p6| for fairly large Ng, 
with which we can compare the VMC results. 

The result of the optimization of variational parameters in t/^ba is shown in Fig. ^ for 
rze = 1/2 as a typical case. The data in the region J/t < 3.3 are fitted to polynomials of 
degree m, where m = 2 for 1/2 and Ug, and m = 3 for ui, respectively. For J/t — * 0, the 
minimum energy is realized for Ui = 1/2 = and ^ 2. In the case of finite J/t, the optimal 
variational parameters show the coupling of charge and spin degrees of freedom as expected, 
i.e., 1^1 7^ and/or 1^2 7^ 0. z/2 and z/^ decrease with J/t while the dependence of z/i on J/t 
is week for J/t < 2. Near J/t = 2, Us intersects the zero line, and z/2 becomes —1. For 
larger J/t, the attractive correlation between electrons with opposite spins is prominent. 
For J/t > 3.3, the variational state separates into the electron rich and poor phases. The 
variational parameters abruptly change their behavior at the phase separation boundary as 
shown in Fig. |I|. The data in the phase separated region are fitted to other polynomials. 

For other values of rie, the optimal values of variational parameters show similar behav- 
iors except for curvatures. Notice that the trial state with optimal variational parameters 
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is singlet although the Jastrow factor Fj is spin-dependent; The VMC evaluation of the 
physical quantities shows that the total spin is zero, and the spin correlation function is 
isotropic {S^^{k) = S^^^k) = S^^{k)), as far as the optimized ipBA is employed. 



Next we discuss the variational energies. It has been shown that the ground-state energy 
converges smoothly to the thermodynamic limit. JT7|JTP|] Following Yokoyama and Ogata, 



l7| , p!9|| we estimate the variational energy per site in the limit Ng ^ oo from the finite size 
scahng. We calculate the variational energies of 12-, 20-, 36-, 60-, and 100- site systems for 
He = 0.5 using the optimized ipBA and ipuu, and then fit the results to the formula 

E/N, = E^ + Cr/Nl + C^/Nt + C,/Nl (18) 

The fitted values of -E'oo(V'ba) and -Eoo('^hm) are listed in Table H for several values of J jt. 
They are compared to the exact results obtained from the Bethe ansatz for J/t = p8[ and 



2, p9| and the extrapolated values of the exact diagonalization of small clusters for J/t = 1 
and 3. The latter is evaluated from fitting the energies of 4-, 8-, 12-, and 16- site clusters to 
Eq. (0). H 

Using V^BA, the ground-state energy per site in the limit A^^ — ^ oo is obtained to be 
Eoo = —2t/iT for J/t ^ 0, equivalent to the exact energy. The reason for the coincidence 
is that the energy is determined only by the charge degree of freedom in the limit J/t ^ 0, 
whose treatment is rigorous in V'ba- In fact, the variance {{H — E{;ip-Qx)Y) in VMC sweeps 
vanishes at any rie- For all the range of J/t, -E'oo(V'ba) is quite close to the exact energy, 
as shown in Table |. Especially in the small J / 1 region, the advantage of ■j/'ba over ipmA is 
obvious. The difference in energy between V^hm and the exact one is largest for J/t = 0, 
while the error of ipBA gradually increases with J/t except for J /t = 2. 

For J/t = 2, both of -Eoo('^ba) and -Eoo(V'hm) are extremely close to the exact energy 
obtained by the Bethe ansatz. In this connection, Yokoyama and Ogata have shown that 



the Gutzwiller function ?/'g is a good trial function for J/t = 2 and all range of ng. ||T7 
In fact, the analytical calculation of the energy using the Gutzwiller wave function shows 
EooijpG) = —0.903092 . . ., |l30|JT7[| which is close to the exact one. Moreover, Kuramoto and 
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Yokoyama [|3T[ have found that ipc is the exact ground state of the t-J model with long- 
range interactions and hoppings satisfying = 2tij oc r^^^. In the supersymmetric case, the 
ground state behaves as almost a free-electron state [jl9| except for the exclusion of double 
occupancy because of the cancellation of hopping and interacting processes, both of which 
have the same weight t = J/2. Therefore, it is natural to expect that variational energies 
converge to the same value if a trial function can recover the "free-electron" nature for 
J/t = 2. The above three variational wave functions have this property correctly although 
the long-range behavior of correlation functions is different. 

In Fig. ^, (a) the momentum distribution function n{k), (b) the spin correlation function 
S{k), (c) the charge correlation function C{k), and (d) the singlet pairing correlation function 
P{k) are shown for J/t = 0, where open circles denote the VMC results of ipBA, and the 
broken lines represent those of ipnu- We evaluate the data for the lattice with 100 sites 
at quarter filling. The exact results obtained from the Bethe ansatz solution for Ns = 52 



are also shown in Fig. with solid lines. p^p!7|,p6| In the TLL, the momentum distribution 



function exhibits power-law singularities at kp and Skp although the latter is very weak. 
Qualitatively, both of ipBA and i^uM reproduce the anomalous power-law behavior inherent 
in TLL as shown in Fig. |^ (a). However, n{k) by ipnu departs appreciably from that by the 
Bethe ansatz while n{k) by ipBA almost coincides with the exact one. For the correlation 
functions, there are also the apparent differences between ipnu and the exact result while 
ipBA is quite close to the exact one, as shown in Figs. (b)-(d). In particular, C{k) calculated 
with ipUM shows different behavior: it has a small cusp at 2kF. On the other hand, both 
of ipBA and the exact result exhibit linear /c-dependence in C{k) reflecting the free spinless 
fermions nature. Thus ipBA recovers the global features of correlation functions correctly in 
contrast to ipnu- 

In Fig. ^ we show the VMC results of the momentum distribution function and the 
correlation functions for finite J/t and ng = |. We have plotted the results of ipBA (open 
symbols) and ipnu (broken lines) for comparison. At J/t = l(open diamonds), the difference 
between ipBA and ■i/'HM is observed, but it is not so large as in the case of J/t = except for 
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the peak of C(2fcir). For J /t = 2 (open squares), the two variational functions give almost the 
same results, in accordance with the agreement of energies. At J /t = 3 (open triangles), it is 
remarkable that P{k) at A; = is strongly enhanced due to the increased effective attraction 
between neighboring electrons with opposite spins as shown in Fig. ^ (d). One can also see 
the enhancement of P{k) near k = tt as J/t increases. The singularity of P{2kp) becomes 
sharp when J /t = 2. 

The correlation functions by z/'ba shown in Fig. |] agree well with the result of the 
Quantum Monte Carlo calculation or of the exact diagonalization in smaller lattices. 



B. Correlation exponent 

Let us now examine the correlation exponents to discuss the long-range behavior of 
correlation functions. In the TLL regime, both the charge and spin excitations are gapless 
and the correlation functions show power-law decay. Following the TLL theory, the leading 
singularities of the distribution function and correlation functions can be written as follows, 
ii 



n{k) ~ 


\k- 


fci? "sgn(A; — kp) 


for k kp , 


(19) 


S{k) ~ 


\k- 


2kF\'^'^ for k ~ 


2kF , 


(20) 


C{k) ~ 


\k- 


2/^^!"-^ for k ~ 


2kF , 


(21) 


P{k) ~ 




for A; ~ . 




(22) 



Logarithmic corrections have been omitted in these formulae. As far as the interaction is 
isotropic in spin space, the critical exponents are described by a dimensionless TLL param- 



eter Kp as follows. 



a = mm 



1, {Kp - l)7(4i^p)J , (23) 

V = Kp + l, (24) 

^i = l/Kp + l. (25) 
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The correlation exponent Kp can now be rather easily calculated from our VMC result of 
P{k) than n{k)^ C{k), or S{k) since the singularities of these quantities become much weaker 
than that of P{k) in some cases. Following Assaad and Wiirtz, |^ we use the following 
procedure to obtain the exponent of P{k). Let us assume a behavior 

In |P(0) - P{ki)\ = -(/i - 1) InNs + a, (26) 

where ki = 2-71 /Ns, and /i and a are the fitting parameters. Fitting the data for various 
lattice sizes to Eq. (p6|), the exponent of P{k) can be evaluated by finite size scaling. Once 
fi is obtained, Kp is given by Eq. (pSj) . 

As an example, we plot the fitting result of P{k) for = ^ in Fig. ^. We have calculated 
P{k) of 12-, 20-, 36-, 60-, and 100- site systems for = ^ using the optimized V'ba, and then 
fit the results to Eq. (^). The linearity of these plots is good. The correlation exponent Kp 
obtained from the slope of these plots and Eq. ( ]25|) is tabulated in Table |I[ For ipnu, one 
can rather analyze the exponent directly. The relation between the variational parameter in 
ipnM and Kp has been derived analytically from the generalized conformal field theoretical 
argument to be 



Kp = . (27) 

Therefore, the optimal value of z/ can be used to determine Kp for V'hm- The results for ipnu 
in Table H are obtained from this equation. The expected values of the exponents are also 
shown in Table |I|: for J/t = and 2, Kp has been exactly determined from the bosonization 
theory ||3^ or the conformal field theory, [^,0 with Kp = 0.5 and 0.85, respectively, and 



the exact diagonalization on a 16-site ring has shown ^ that Kp ^ 0.62 and 2.2 for J/t = 1 
and 3, respectively. 

The results of ipuM and ipBA are consistent with the expected values; At J/t = 3, Kp 
becomes larger than 1 and thus the superconducting correlations correspondingly dominate 
the long-range behavior while for smaller value of J/t, Kp < 1. However, the quantitative 
coincidence is not so good for ipuM, while the exponents evaluated from ip^A are surprisingly 
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close to the expected values. Therefore, one could conclude that our variational wave func- 
tion ^/'ba can quantitatively reproduce not only the amplitude of correlation functions, but 
also the correlation exponent. The last column in Table |I| shows the value of {ui + 1^2 + 2)~^ 
evaluated from the optimized variational parameters in ipBA- It seems that Kp agrees with 
{ui + z/2 + 2)^^. This point will be discussed in the next subsection. 



C. Phase diagram 

The phase diagram of the ID t-J model obtained by ?/'ba and ^rvb is shown in Fig. ^ with 
contour lines for several values of Kp. As seen in Fig. ^, there are four distinct phases. For 
small J/t, the ground state is a repulsive TLL with Kp < 1. In this region, spin correlations 
dominate the long-range behavior. Increasing J/t, these correlations are suppressed, and 
the ground state changes to an attractive TLL with Kp > 1. It has dominant singlet 
pairing correlations. For larger J/ 1, the variational state is phase-separated into the electron- 
rich phase with antiferromagnetic order and the empty phase. In this region, the longest 
wavelength charge correlation C{ki = 2tt/Ns) diverges when the system size A^^ is increased. 
This behavior is in contrast with the TLL, where C{ki) remains finite. The phase separation 
boundary in Fig. ^ is determined by the behavior of C{ki). These three phases are described 
by "ipBA- The spin-gap state lies between the TLL and phase separated regions at small 
densities, where V'rvb is more stable than tpBA- The optimal values of variational parameters 
in i/jRYB fit to the formulae as h = 1.037-0.151(J/t) and A = 0.781 -0.173(J/t) for Ue = 0.2 
and 2.8 < J/t < 3.2, for example. 

The phase diagram determined by the present method is consistent with the result of 
the exact diagonalization on a 16-site ring. p[ It is correctly recovered that Kp = 1/2 at 



any electron density in the limit J/t 0. p^,|35[] When we compare the spin-gap region 



with the result of the wave function of noninteracting singlet pairs by Chen and Lee, 
the correlation effect of the Jastrow factor introduced in V'rvb seems to play no essential 
role since these two trial functions give almost the same result; The correlation effect only 
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slightly pushes the spin-gap region to larger values of rie. This is because the balance of 
energy hardly changes, i.e., the magnitude of correlation energy induced by the Jastrow 
factor in ipBVB is less than 2%, which is comparable to the energy lowering of tpBA compared 
to ipuu- However, there is a significant difference in the long-range behavior of correlation 
functions. The wave function by Chen and Lee gives Kp = oo while Kp is finite for our ipRVB 
as seen below. 

Figure ^ shows the distribution function and correlation functions for 150 sites and 30 
electrons. The selected values of J/t are 0, 2.5, and 3 as typical cases of the repulsive 
TLL, attractive TLL, and spin-gap state, respectively. Triangles and squares in Fig. ^ are 
evaluated with ip^A while circles with V'rvb- The spin correlation functions for J/t = and 
2.5 exhibit the linear behavior at small k characteristic of TLL as shown in Fig. ^ (b). On the 
other hand, S{k) for J/t = 3 is quadratic at small k and analytic for all wave vectors. Unlike 
TLLs, a Luther-Emery liquid exhibits exponential decay of the spin correlation function in 
real space, while both charge and singlet pairing correlations decay algebraically. ^ The 
short-range behavior of S{k) for J/t = 3 shows that the spin-gap state is characterized 
as a Luther- Emery liquid. The charge correlation function for J/t = 3 exhibits a cusp at 
k = 2k p, indicating the formation of bound singlet pairs. As J/t decreases, this cusp is 
suppressed. More definitive characters of the three phases can be seen in the singlet pairing 
correlation functions plotted in Fig. ^ (d). P{k) is fully suppressed when J/t = 0. As J/t 
increases, P{k = 0) becomes much larger, indicating the growth of long-range order. The 
cusp at /c = is greatly enhanced for J/t = 3. 

Finally, we plot the exponents of the singlet pairing correlation function for He = 0.2 as 
a function of J/t in Fig. |^, where /x — 1 is evaluated from fitting P{k) to Eq. ( p6D for lattice 
sizes ranging from A^^ = 30 to 150. In the region 2.8 < J/t < 3.2, the variational energy of 
^RVB is lower than that of ^ba- Correspondingly, squares and circles in Fig. |^ are evaluated 
with ipBA and ipm/B, respectively. In addition, z^i -|- z/2 + 2 for ipBA and 2A — 1 for ipByB are 
shown in Fig. |^ with solid lines, where z/i, i>2, and A are optimized variational parameters. 
It seems that the fits of z/i + z/2 + 2 and 2A — 1 to the exponents of P{k) are good. The 
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exponents of the correlation functions in the Luther-Emery hquid that decay with power 
laws can be also described by a single parameter Kp like TLL, ^ and it holds that 

^i = ^. (28) 

Using Eqs. (^) and one can conclude from Fig. |^ that Kp relates to the variational 
parameters in -j/^ba and V'rvb as 

Kp = ^_ for ^BA, (29) 

Z/l + Z/2 + 2 

and 

= ^ for ^RVB, (30) 

respectively. The critical exponent Kp evaluated from P{k) for Ue = 0.2 is tabulated in 
Table |in| together with the predicted values by Eqs. and (ppp, and they agree very 
well. For other values of rig, Kp also agrees with {ui + 1^2 + 2)^^ or (2A)~^. The example for 
rie = 0.5 is shown in Table Kp shown in Fig. |^ is evaluated from Eq. (|29|) . 

The relations (pOf ) and (^Oj) are confirmed by the following facts, (i) In the limit ne 1, 
we have z^i = z/2 = and Us ~ 2, i.e., ipBA = Y for the Heisenberg chain since X becomes 
only a constant. This is in accordance with Kp 1/2 as Ue —* 1. (ii) It is correctly 



recovered that Kp 1/2 at any electron density in the limit J/t 0, p^j35[1 where we 
have Ui = 1^2 = 0. (iii) The magnitude of the discontinuous jump in /i — 1 at J/t ~ 2.8 
is (z/i + U2 + 2) — (2A — 1) ^ 1 as seen in Fig. |^ This is consistent with the crossover 
from the TLL to the Luther-Emery liquid behavior described by Eqs. ( pQ]) and (|30D. (iv) 
^1 + ^2 + 2 ~ on the phase separation boundary in Fig. ^, which leads to Kp = 00. This 
corresponds to the divergence of C{ki). 



IV. SUMMARY AND DISCUSSION 



In this paper, we have carried out the VMC calculation for the ID t-J model. As a trial 
state, we have proposed a new type of variational wave function based on the Bethe ansatz 
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solution. In this wave function, the separation of charge and spin degrees of freedom is real- 
ized explicitly, and the long-range correlation factor of Jastrow-type is included. With this 
wave function, it has been shown that the remarkable improvement is achieved especially in 
the small J/t region: the variational energy, momentum distribution function, and various 
correlation functions exhibit an excellent coincidence with exact ones. The evaluation of 
correlation exponents with the finite size scaling has shown that this variational wave func- 
tion can correctly reproduce not only the global features of correlation functions but also 
the long-range behavior with anomalous power-law decay, which is characteristic of TLL. 

In addition, a variational wave function of singlet pairs correlated with a Jastrow factor 
has been introduced to describe the spin-gap phase. This wave function correctly exhibits 
enhanced superconducting correlations and exponential decay of the spin correlation func- 
tion, as expected for the generalized Luther-Emery state. 

Comparing the energies of the trial function based on the Bethe ansatz solution and the 
generalized Luther-Emery state, the whole phase diagram has been determined. The VMC 
results show that our wave functions provide a more precise description of the ground-state 
properties for the ID t-J model in the whole phase diagram. Evaluating the correlation 
exponents by the finite size scaling analysis, the relations between the exponent Kp and the 
variational parameters in the trial functions have been established. 

Let us now compare our trial wave functions with others. For strongly correlated elec- 
tron systems, the Gutzwiller- Jastrow-type trial state has been extensively studied, but the 
conventional Jastrow factor does not recover the expected TLL behavior if only short-range 



correlations are included. |ll7| , p!9| A trial wave function of this type is essentially a Fermi 



liquid state. The wave function introduced by Hellberg and Mele successfully exhibits the 



power-law singularity of TLL. |T8| The long-range nature of the Jastrow factor is essential 
for the non- Fermi liquid behavior. However, the correlation exponent does not coincide with 
the exact value. This disagreement becomes apparent when we compare the global features 
of various correlation functions. The difference between ipmn and the exact result is largest 
for J/t = 0. These are in sharp contrast to ipBA, with which we can quantitatively reproduce 



both of the amphtude and the exponent of correlation functions. This is because i/jba has 
the separation of charge and spin degrees of freedom correctly. The effect of the spin-charge 
separation becomes clear especially in the small J/t region. In fact, when we compare the 
phase diagram determined by the present method with that of ipnu, [0 the behavior of Kp 
is much improved in the repulsive TLL region, while it agrees with the result of ipuu in the 
attractive TLL phase. 

Finally we mention some remaining issues. The spin wave function Y{yi, . . . ,yM) in 
Eq. (1^) has been approximated as a trial function of Jastrow-Marshall-type for its simplicity. 
It is known, however, a liquid state is realized in the ID antiferromagnetic Heisenberg model. 
When a RVB-type trial state is used as the spin part in Eq. (|^ instead of Jastrow-Marshall- 



type, further improvement may be expected. Actually, Ogata has examined a trial 
function composed of a spin trial state of RVB-type and spinless fermions for J/t = 0. It 
is interesting to correlate this trial function with a Jastrow factor, applying to all range of 
J/t. 

An application of present method with some modifications to magnetic properties such 
as spin susceptibility and magnetization curve is also interesting in order to elucidate the 
metal-insulator transition in the ID t-J model. In fact, it was shown that Jastrow wave 
functions reproduce charge and spin susceptibilities and magnetization curve correctly, in 
contrast with the Gutzwiller approximation. [|19| 

An important question is whether the properties of strongly correlated electrons realized 
in ID system can be extended to higher dimensions because of their close connection to 
high-Tc superconductors. |jl2[ In 2D system, it is not established even for the metallic regime 
whether the ground state is the Fermi liquid or TLL. In these contexts, an extension of the 
present method to 2D systems together with reexamination of the appropriate Hamiltonian 
is under consideration. 

Quite recently, we have found that a part of our result for the variational wave function 
of correlated singlet pairs has been obtained independently by Chen and Lee. 
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FIGURES 

FIG. 1. Optimization result of variational parameters in ipBA is shown for Ug = ^ and 
Ns = 100. Triangles, squares, and circles represent z^i, 1/2, and z^^, respectively. The solid lines are 
the least-squares fits of the data. The transition to a phase-separated state is shown by an arrow 
at J/t w 3.3. 

FIG. 2. (a) the momentum distribution function n(k), (b) the spin correlation function S{k), 
(c) the charge correlation function C{k), and (d) the singlet pairing correlation function P{k) for 
Ue = and J/t = 0. Open circles and broken lines denote the VMC results of iI^ba and VhMj 
respectively, for Ns = 100. The exact results 

FIG. 3. (a) the momentum distribution function n{k), (b) the spin correlation function S{k), 
(c) the charge correlation function C{k), and (d) the singlet pairing correlation function P{k) for 
Ue = and Ns = 100. The data are evaluated by the VMC calculations with ipBA for J/t = 1 (open 
diamonds), 2 (open squares), and 3 (open triangles). The broken lines denote the corresponding 
results of tpuM- 

FIG. 4. The VMC results of In |P(0) - P{ki)\ at rie = 1/2 are plotted for J/t = (circles), 1 
(diamonds), 2 (squares), and 3 (triangles) as a function of InA'^s, where ki = 2tt/Ns. The data are 
evaluated with ^ba for Ns = 12, 20, 36, 60, and 100. The solid lines are the least-squares fits of 
the plots. 

FIG. 5. The phase diagram of the ID t-J model as determined by Vba and V'RVB- In the 
spin-gap phase, ^/'rvb is more stable than ^ba while other phases are described by iI^ba- The curves 
represent the contours of constant correlation exponent Kp evaluated from Eq. (]29|). 

FIG. 6. (a) the momentum distribution function n{k), (b) the spin correlation function S{k), 
(c) the charge correlation function C{k), and (d) the singlet pairing correlation function P{k) for 
He = 0.2 and Ns = 150. Selected values of J/t are (triangles), 2.5 (squares), and 3 (circles) as 
typical cases of the repulsive TLL, attractive TLL, and spin-gap state, respectively. 



10,17,2q] for Ns = 52 are also shown with solid lines. 
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FIG. 7. Exponents of the singlet pairing correlation function for ?ie = 0.2 as a function oi J/t. 
Squares are evaluated with V'BA) while circles with i/'rvBi fitting P{k) to Eq. (|^) for lattice sizes 
ranging from Ng = 30 to 150. Solid lines represent the value vi + U2 + 2 for -i/^BA and 2A — 1 for 
V'RVB, where i^i, 1^25 and A are the optimized variational parameters. 
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TABLES 

TABLE L Ground-state energies of the ID t-J model in the hmit Ng ^ oo for the quarter-filled 
case (ne = i). The VMC results are obtained from the finite size scaling Eq. (|l8| ) with Ng = 12, 
20, 36, 60, and 100. Exact results obtained from the Bethe ansatz for J = and 2, and 
extrapolated values of the exact diagonalization of small clusters for J = 1 and 3 are listed as 
a comparison. The unit of the energy is t. 



J 



Exact 



VMC ^oo(V'hm) 



(Error 



VMC ^oo(V'ba) 



(Error %] 



-0.636620 -0.6119 ±0.0003 (4) 

-0.755359 -0.7493 ±0.0001 (0.8) 

-0.903649 -0.90315 ± 0.00005 (0.06) 

-1.081713 -1.0774 ±0.0001 (0.4) 



-0.636620 ± (0) 

-0.75488 ±0.00005 (0.06) 

-0.90318 ±0.00006 (0.05) 

-1.0806 ±0.0001 (0.1) 



TABLE IL The correlation exponent Kp for Ue = 2- Results for ipBA are evaluated from fitting 
P{k) to Eq. ( p^ ) for lattice sizes ranging from A''^ = 12 to 100. With V'hMi Kp is directly obtained 
from Eq. (p7[). The last column is evaluated from the optimized variational parameters in tpBA- 





Expected 








J/t 


value 


Kp{il)mA) 


^p(^ba) 


(i/i ± ± 2)-i 





0.50 |34,3| 


0.40 


0.51 ±0.03 


0.50 


1 


0.62 i 


0.57 


0.62 ±0.05 


0.65 


2 


0.85 [36| 


0.95 


0.85 ±0.07 


0.95 


3 


2.2 II 


3.2 


2.2 ±0.2 


2.4 
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TABLE III. The correlation exponent Kp for Ue = 0.2, obtained by fitting P{k) to Eq. (|26|) 
for lattice sizes ranging from Ng = 30 to 150. The last column is evaluated from the optimized 
variational parameters in the trial wave functions. 



(i/i + V2 + 2)-^ 



J/t Type of Kp or {2\)-^ 

0.0 BA 0.51 ±0.05 0.50 

1.0 BA 0.57 ±0.07 0.63 

2.0 BA 0.97 ±0.06 0.98 

2.5 BA 1.7 ±0.2 1.6 

2.8 RVB 1.6 ±0.2 1.7 

3.0 RVB 1.8 ±0.2 1.9 

3.2 RVB 2.3 ±0.4 2.2 
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